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ABSTRACT 

A  hybrid  mode  analysis  of  microstrip  on  a  dielectric 
substrate  is  presented.   A  numerical  solution  for  wavelength 
and  characteristic  impedance  of  single  and  coupled,  balanced 
strips  is  obtained.   Line  parameters  are  shown  to  be  very 
frequency  dependent  but  in  agreement  with  the  quasi-static 
results  of  other  investigators  in  the  low  frequency  limit. 

The  extension  of  this  technique  to  perturbation  analysis 
of  microstrip  on  ferrite  is  also  described.   A  numerical 
solution  for  the  propagation  constant  is  obtained. 
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I.   INTRODUCTION 

The  analysis  of  micro-strip  is  of  great  importance  as 
this  type  of  transmission  line  has  found  wide  use  due  to 
its  compatibility  with  microwave  integrated  circuitry.   The 
open  microstrip  transmission  line  has  been  predominant  be- 
cause it  can  be  etched  or  deposited  easily  on  the  substrate, 

The  early  analytic  work  on  microstrip  lines  has  been 
based  on  a  proposed  TEM  mode  of  propagation  which  is,  in 
essense,  a  static  approximation  to  a  dynamic  system.   For 
sufficiently  low  frequencies  the  quasi-static  theory  can 
be  employed  to  obtain  the  characteristics  of  microstrip 
lines.   When  the  wavelength  in  a  microstrip  line  becomes 
comparable  to  transverse  dimensions  of  the  line  the  devia- 
tion from  quasi-static  behavior  becomes  significant  and 
high  order  modes  of  propagation  become  possible. 

In  this  present  work  a  Spectral  Domain  transform  method 
is  applied  for  calculating  the  frequency  dependent  charac- 
teristics of  single  and  coupled  microstrip  transmission 
lines.   Effects  of  geometry  on  the  dispersion  and  charac- 
teristic impedance  have  been  analyzed. 

Figures  1  and  2  show  the  geometry  of  the  microstrip 
line.   The  strip  conductor  is  assumed  to  be  infinitely 
thin  and  perfectly  conducting  and  the  ground  plane  and 
dielectric  or  ferrite  substrate  are  assumed  to  be  infinite 
in  extent. 
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Figure  1.   Three  Dimensional  Microstrip  Transmission  Line 


11 


igzzzsa {zzzzzza 

r  W  -*^-   S  "+■  W  -*• 


Region    1 

w  //  /  z  //i  KSZZZZZ1 


Figure  2.   Coupled  Microstrip  and  Regions 
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It  is  assumed  initially  that  the  substrate  material  is 
lossless  and  isotropic.   Lossy  anisotropic  substrates  are 
treated  later  using  perturbation  theory. 

The  Spectral  Domain  transform  method  was  suggested  by 
Itoh  and  Mittra,  [Ref.  6]  and  yields  a  solution  to  the 
Boundary  value  problem  of  microstrip  via  the  method  of 
moments.   The  method  was  used  by  Itoh  and  Mittra  to  find 
the  dispersion  characteristics  of  a  single  strip.   In  this 
work  the  method  has  been  extended  to  cover  dispersion  and 
characteristic  impedance  of  both  single  and  coupled  micro- 
strips.   The  case  of  a  microstrip  on  a  ferrite  substrate 
has  also  been  solved  using  perturbation  theory. 
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II.   DISPERSION  CHARACTERISTICS  ON  DIELECTRIC  SUBSTRATE 

A.   FIELD  AND  BOUNDARY  CONDITIONS 

Let  the  following  electric  and  magnetic  fields  exist  and 
propagate  in  the  z-direction  as  justified  by  the  Hertzian 
Vector  Potential  functions  described  in  Appendix  A. 


H   = 


1  2aG  YZ 

k  <b    e  ' 

cY 


i  2  j.e  YZ 
k.  <b  e  ' 
c 


(1) 
(2) 


where  y    is  the  Propagation  constant.   For  the  lossless  sub 
strate  material 

Y  =  j3. 


From  the  field  expressions  in  equation  (1)  and  (2)  all  other 
components  of  electric  and  magnetic  field  can  be  derived 
from  Maxwell's  curl  equations,  as  shown  in  Appendix  B  and 
leads  to  the  following  equations. 


H. 


H 


Y 


Ml 
9x 

h 


jcou 


9c}> 


h\ 


9y 


9cf>       J      .  3<|> 

,       9x         J  9y 


.3* 


9$ 


.YZ 


.YZ 


.YZ 


9y         J  9x 


(3) 

(4) 
(5) 
(6) 


1  4 


Applying  boundary  conditions  at  the  interface  between 
region  2  and  the  ground  plane,  tangential  electric  fields 
must  be  zero,  so  it  follows  that  at  y  =  0, 


E   (x,0,z)  =  0 
z2 


(7) 


E   (x,0,z)  =  0. 
x2 


(8) 


Also  at  the  interface  between  regions  1  and  2,  tangen- 
tial electric  fields  must  be  continuous: 
At  y  =  D, 

(9) 


E   (x,D,z)  =  E   (x,D,z) 
zl  z2 


E   (x,D,z)  =  E   (x,D,z) . 
xl  x2 


(10) 


Also  the  electric  fields  will  exist  only  in  the  dielec 
trie  part  of  the  interface  and  can  be  expressed  as, 


0 


on  strip 


E   (x,D,z)  =  { 


Y  Z 

e    (x) e '    elsewhere 

z  v  J 


(11) 


E   (x,D,z)  =  I 
Xl 


on  strip 


Y  Z 

e    (x)e'    elsewhere 

x  *-  J 


(12) 


Similarly,  tangential  magnetic  fields  must  be  discon- 
tinuous by  corresponding  surface  current  densities 


H   (x,D,z)  -  H7  (x,D,z)  =  ( 
zl  ■  ^2 


Y  Z 

J  (x)e'   on  strip 


elsewhere . 


(13) 
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'  J    (x)  eYZ 


H   (x,D,z)  -  H   (x,D,z)  =< 


0 


on  strip 


elsewhere 


(14) 


Substituting  the  field  expressions  of  equation  (1) 
through  (6)  into  the  boundary  condition  expressions  of 
equations  (7)  through  (14)  ,  one  obtains 


c2  L 


=  0 


(15) 


y~tt    (x,0)  -  ja3y2  -^ 


(x,0)  =  0 


(16) 


k*  <fr®(x,D)  =  k*  4>e(x,D) 
cl  c2 


(17) 


Y- 


3^(x,D) 


3x 


jwy 


3 <()!,_  (x,D) 

5y 


=  Y- 


34>|(x,D) 


3x 


jo)y. 


34>2(x,D) 


(18) 


k*  *l(x,D) 
cl  -1 


e    (x) 


on  strip 


elsewhere 


(19) 


Y- 


34>1(x,D) 


3x 


jtoy 


1  ""  9y 


k*  <£(x,D) 
cl 


k*  <^(x,D) 


on  strip 


e  (x)  elsewhere 

x^  ' 

J  (x)  on  strip 


elsewhere 


(20) 


(21) 
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Y- 


9<f>1(x,D) 


3x 


34>2(x,D) 

9x 


+  Jwei  — §7— } 


3<J>|(x,D) 

+  Jwe2  — 3y 1 


J  (x)   on  strip 


(22) 


0      elsewhere 


B.   SPECTRAL  DOMAIN  TRANSFORM 

All  potential  functions  must  satisfy  the  following  rela 
tion: 


V24>  +  k2  4>  =  0 
t     c 


(23) 


where 


k2    =    2   +  k2   =  k2   _   32 


(24) 


One  introduces  the  Fourier  transform  Lo  the  a -domain,  as 
suggested  by  Itoh  and  Mittra: 


+  oo 


1  ax 


*i(a,y)  ■  j       4)i(x,y)eJV*"dx    i  =  1,2,3...     (25) 


Thus,  the  transform  of  equation  (23)  becomes 


^x  i  0  i  +  Fx  t  0  ]  *  K  ^x  [  *i  ■  o 


(26) 


(ja)2  F   [  <j>(x,y)  ]  +   «   F   [  <j)(x,y)  ] 


By2   x 


+  k2  Fx  [  c|>  ]  =  0 


where 


x  L 


9*(x,v) 


9x 


1  = 


-ja  Fx    [    <J)(x,y)  ] 


has  been  applied. 


(27) 


(28) 
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Clearly,  one  can  obtain  a-domain  representation  of 
equation  (23) , 

-a2<Ka,y)  +  -£t   *(<*»y)  +  K      *(<*.y)  =  0 


dy 


(29) 


or 


iigSUTl-  (a»  -  ki)  t(a,y).  (30) 

Equation  (30)  should  be  analyzed  carefully  both  for 
Region  1  and  Region  2. 

For  Region  1,  (free  space  propagation  constant), 


•  2_^2 


where 


a2  -  k2   =  a2  +  62  -  k! 


jo     271 

and  3  =  t~t 


(31) 


k1  =  co/y  e 
1      o  o 


2       2      /-   2  TT   -v  9 

YJ_  =  a2  +  (  ~  )2 


t^)2-  i 


(32) 


X'  is  the  microstrip  wavelength  and  is  related  to  the  free 
space  wavelength  X  as 


X  >  X'  > 


(33) 


By  substituting  Equation  (33)  into  (32)  one  can  find 
lower  and  upper  values  for  y-i 


a2  +  (er   -  1)(  2£  )2  >  y^  >  a2. 


(34) 


Therefore  Yi  is  always  a  real  quantity  independent  of 
the  values  of  a. 
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Similarly  for  Region  2, 


2       2 


C2J,. 


(fr) 


(35) 


so  it  is  clear  that  y?   will  be  imaginary  for 


27T 


2tt 


,er  "  (  >j   )2  <  a  < 


-  (  t4)2 


2  A2 

and  will  be  real  when  a  is  in  the  range 


■°°<a<  - 


2tt 


T\/er, 


V.   "\  I   J        j 


X 


X 


X 


(36) 


(37) 


(  jj   )2  <a<+~. 


For  Region  1,  the  field  function  differential  equation 
solution  is  as  follows: 


8y2   1 


(a,y)  =  Y2(Ka,y)  =  (az-k*  )$1(a,y) 


Cl   X 


(s2  -  Y2)$i(s)  =  0. 


1j^1 


(38) 
(39) 


Therefore , 


0. 


(a,y)  =  R1(a)e  +  Yiy  +  R2(a)e"Yly. 


(40) 


Clearly,  for  infinite  positive  values  of  y,  the  field 
should  vanish,  or, 


lim  <f>  (a,y) 

y-voo 


=  0 


(41) 


then, 


x(a,y)  =  R2(a)e"Yly  =  A(a)e"Yly 


(42) 


where 
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A(a)  ■  R2(a). 

Similarly  for  Region  2,  there  are  two  solutions,  cor- 
responding to  the  real  or  imaginary  character  of  y?.  The 
solution  for  y?  imaginary  is 


$2(a>y)  =  B(a)SinY2y  +  C(a)CosY2y 


(43) 


where 


Y?  =  JY! 


and  for  y~  real  it  is 


$2(a>y)  =  B(a)Sinhy?y  +  C(a)Coshy~y 


(44) 


e        h        e        h 
Therefore,    the    transforms    of   <j>- ,    <{>..,    <J>?,    <j)?    are    given 


by 


Region  1 
*®(a,y)  -  Ae(a)e 

h  h       -Yi(y-D) 

$*(a,y)  =  An(a)e   1 


(45) 
(46) 


Region  2 


B^(a)SinhY2y  +  C®(a)Coshy2y  ,  Y2  Real   (47) 


$2(a,y)  =  < 


(48) 


jB„  (a)Siny!'y  +  CT(a)Cosy"y   ,  Y?  Imaginary 


BH(a)Sinhy2y  +  C^(a)Coshy2y  ,  Y2  Real   (49) 


4>?(a,y)    -   < 


h  h  C50) 

jBT(a)SmY2y    +    (^(a)  Cosy'^y       ,    Y2    Imaginary 
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where  superscript  (e)  denotes  the  electric  field  case  and 
(h)  magnetic  field  case  and  where 

^   '  Y2  =  "JY2- 
Spectral  Domain  representation  of  the  boundary  condition 
expressions  can  be  obtained  by  taking  the  Fourier  transform 
of  the  equations  (7)  through  (14)  with  the  following  results: 

k>  *f (a,0)  =  0  l  (51)  ' 

-jaY*®(a,0)  "  >y2  -g±  (a,0)  =  0  (52)  ' 

k*  $e(a,D)  =  k2    *?(a,D)  (53) 

cl  c2  * 

8^(a,D) 
-jaY*1(a,D)  -  ja>y1 =  -jay$2(a,D) 

"  >y2  97 

k*  $®(a,D)  =  *z(c0  (55)  • 

-jay<^(a,D)  -  jcoyi  -^  (a,D)  =  ^(a)  (56) 

k*.*!j(cx,D)  -  k2c    »5(afD)  ■  ^(a)  (57)  ^ 


jay$1(a,D)  +  jwe1  ^- +  jaY$2(a,D) 


e 


^(a,D)  (58) 

■>£2  87 =  JzM' 

Where  the  derivative  transform  pair 

F    j  UfeXL   ]  .  -ja#(a.y) 
has  been  applied. 
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By  substituting  the  Field  expressions  for  real  (hyperbol- 
ic case) and  the  imaginary  (trigonometric  case)  values  of 
Y?  into  the  transformed  equations  (51)  through  (58) ,  the 
following  equations  are  obtained. 

Hyperbolic  Case: 

k*  C®(a)  -  0  .  (59) 

c2  n 

-j[ayC^(a)  +  a3y2Y2B^(a)]  =  0  (60) 


k*    Ae(a)    =    k*     [C®(a)CoshY2D    +    B®  (a)  Sinhy^]  (61) 


V 


j(ayAe(a)    -    toy^A    (a))    =    -    j  ay  [B^(a)  Sinhy^    + 


(62) 


CH  (a)Coshy7D]     -    ju)y?y9  [B„  (a)  Coshy9D    +    C„  (a)  Sinhy^D] 
k*    Ae(a) "-    Bz(o)  (63)     / 

-j(ayAe(a)    -    ajy^A^a))    =    ^(a)  (64)     / 

k*    Ah(a)    -    k*     [Bjj(a)SinhY2D    +    cj]  (a)  Coshy-^D]    =    Jx(a)(65) 

-j(ayAh(a)    +    coe^A^a))    +    j  ay  [Bjj  (a)  Sinhy2D    +    cjj  (a)  Coshy2D] 

-ja)e2Y2[BIeI(a)CoshY2D    +    c|j  (a)  Sinhy^]    =    J    (a)  .  (66) 

Similarly    for    imaginary   values    of   y~   by  making    the    sub- 
stitutions 

Y2    =    JY'2' 
Sinhy2D   =    jSiny"D  (6  7) 
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CoshYoD  =  Cosyp 

in  equations  (59)  through  (66)  a  corresponding  set  for  the 
trigonometric  case  can  be  obtained. 

Finally  constants  for  electric  and  magnetic  field  expres 
sions  are  determined  as  follows: 


Ae(a)  =  pi-  2?  (a) 
Cl 


(68) 


Ah(a) 


ay 


jcoy1Y1   x 


E    (a)  +  "TT  ,  E    (a) 


wylYl~c. 


(69) 


CT  H^  =  ° 


BT,H(C°  =  ° 


e  1 

Brp(a)  =  ^t-7 — c~- nrr  E    (a) 

Ty    J         jk2  Siny'-'D  zy    J 

c2 


Bu(°0  =  i-? — c-= — u n  E     (°0 

H v  J         k2  Sinhy-D  zK    J 


(70) 
(71) 

(72) 
(73) 


CH(a)  = 


ooy?YoSinhY?D 


ay 


E    (a) 

z  v    J 


j^x(a) 


(74) 


CT(a)    = 


oovUYTSinYo'D 


ay 


B  fa)  -  j£  (a) 


k2    zv"'    ■*  x 


(75) 


The  subscript  T  denotes  the  trigonometric  case  and  H 
denotes  the  Hyperbolic  case.   Substituting  equations  (70) 
through  (75)  into  equation  (65)  and  (66)  one  can  obtain  two 
sets  of  equations  of  the  form 
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F1(a,3)£'x(a)  +  F2(a,3)  B2(o)  =  Jx(a) 


F3(a,3)ffx(a)  +  F4(a,3)  i?z  (a)  -  «7x(a) 


(76) 
(77) 


where 


Y  =  j3. 


For  the  hyperbolic  case  the  F.  have  the  form 


F,     (a, 8)    = 


k2 
c 


k2 
c. 


ojyy-L        00]J2Y2 


Cotanhy-D 


(78) 


F   (a, 8)  = 


a8 


wy^ 


— ^£-  Cotanhy9D 

o)y2Y2        r2 


(79) 


F,  («,B)  =  -  F   (a, 8) 


f801 


F4  (a, 6)  =  -  j 


,  Ds2    wej-,    CotanhYoD 

C2 


03y  1  Yi  kT 

11  C-.     c. 


(a3)2 

- — — —  -  coe0Yo 

wy2Y2      2' 2 


(81) 


A  similar  set  of  equations  for  the  trigonometric  case 
can  be  obtained  by  use  of  equation  (67) .   Solutions  of 
equations  (76)  and  (77)  for  E    (a)  and  E    (a)  in  terms  of 
J    (a)  and  J    (a)  may  be  obtained  as  follows : 


EXW    - 


F4(a,8)Jx(a)  -  F2(a,B)Jz(a) 


J) 


N 


(82) 
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-F-(a,BHv(cO  +  F  (a,3V  (a) 
*  (a)  =  — *- g i 5 (83) 

where 

DN  =  F1(a,3)F4(a,3)  -  F2 (a, 3) F3 (a, 3) .  (84) 

Define  the  following  constants, 

F,(a,B) 
GiCo.B)  =  -5-u (85) 

-F  (a,3) 
G?(a,3)  =  — ^ (86) 


F3(a,B) 


G3(a,3)  =     p  (87) 


G4(a,3)  =  -^ (88) 


'N 
finally  obtaining 


Fx(a)  =  G1(a,3)e/X(a)    +  G2  (a,  3)<7Z  (a)  (89) 


2?z(a)  =  G3(a,3)Jx(a)    +  G4 (a, 3) Jz  (a) .  (90) 


C.   PHYSICAL  PARAMETERS 

General  equations  (89)  and  (90)  show  no  dependence  on 
the  actual  physical  configuration  other  than  the  boundary 
conditions . 

In  general,  the  transformed  surface  current  densities 
J"  (a)  and  Jn  (a)  can  be  expanded  in  a  set  of  known  basis 
functions  as, 


J    (a)  =  .E-  a.  j  .  (a)  (91) 

xv  '    1=1   iJxiv  '  *■  J 

J    (a)  =  .if-  b.  j  .  (a)  .  (92) 

z v  J         1=1   i J  zi v  '  v   ' 

Since  the  two  conductors  are  sufficiently  narrow,  one 
may  assume  that  the  surface  current  density  in  the  x-direc- 
tion  is  zero.   By  modifying  equations  (89)  and  (90),  also 
substituting  (91)  and  (92)  one  obtains 

G2(a,3)  i|1  b.jzi(c0  =  *x(a)  (93) 

G4(a,B)  jl1   bijzi(a)  =  »z(o).  (94) 

Furthermore,  although  one  could  perform  a  multiterm  ap- 
proximation to  the  surface  current  density  in  the  z-direc- 
tion,  in  this  case  a  one  term  approximation  was  found 
sufficient  after  different  distributions  were  tried.   For 
the  numerical  results  of  this  present  work  uniform  surface 
current  density  is  used,  thus 

G2(a,6)b1jz.  (a)  =  Fx(a)  (95) 

G4(a,B)b1Jz. (a)  =  ffx(a) ,  (96) 

or  substituting  equation  (92)  for  the  case  where  i  =  1, 

G2(a,3)Jz(a)  =  2?x(a)  (97) 

G4(a,3)Jz(a)  =  *z(a).  (98) 
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It  is  clear  that  one  only  needs  to  work  with  either  one  of 
the  above  equations.   Assuming  a  uniform  current  distribution 


J  (x) 

z  v  J 


-< 


I  /W  on  the  strip 


elsewhere 


where  W  is  the  width  of  the  strip  and  I   is  the  magnitude 

r      o  6 

of  the  current,  one  can  obtain  the  a-Domain  Fourier  trans- 
form, J    (a) .   This  was  done  for  three  different  configura- 
tions, namely  single  strip,  coupled  strips  even  mode  and 
coupled  strips  odd  mode. 

For  the  geometry  in  Figure  3  one  obtains  the  following 

Single  Strip: 


W/2 


J    (a) 
z  v  J 


0_ 

W 


W/2 


eJaxdx,  «Tz(cO 


=  I 


SinC-j) 

aW 

2 


(99) 


Coupled  Strips,  Even  Mode 


:       f-S/2  (S/2+W) 

JZW    =  -r-     J      eJ   dx  +   /      eJ   dx 
-(S/2+W)  ^S/2 


j  (a)  .  2i  cos  SlflWI 
z  v  ^      o      2 


Sin(  -j-   ) 


(100) 


Odd  Mode 


V°)  =  "¥ 


.-s/ 


(S/2MV) 


j  ax 


dx 


(S/2+W) 


S/2 


jax 


dx 


(101) 
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t  J  (x) 


0_ 

W 


■**  X 


W 
2 


W 
2 


*J  (x) 

z  ^  J 


ta  I 
T  o 


W 


■*-  x 


S-W   "S 


2    2^ 


A  J  (x) 

zK    J 


_0 

W 


S/2 


|+W 


>•  x 


2  n 


I  +■ 
o 

W 


Figure  3.   Single  and  Coupled  Strips  liven  and  Odd  Longitudinal 
Current  Distributions. 
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ars+Wl         Sin  (ay) 
J    (a)    =   J2I   Sin  2Li^JLL       ^L_   .  (10i) 

1      2    J 


It  is  convenient  to  make  equation  (98)  independent  of  E    (a) . 
This  can  be  done  by  taking  an  inner  product  with  a  function 
orthogonal  to  E    (a).   Applying  this  concept,  equation  (98) 
becomes 


<G4(a,eHz(a),W(a)>  =  <2?z  (a)  ,W(a)  >  (102) 

where  the  inner  product  is  defined  by 

<Y(a),W(a)>  =    f        Y(a)W(a)da.  (103) 

-co 

A  suitable  weighting  function  W(a)  is  the  complex  con- 
jugate of  J    (a);  i.e.  J*  (a)  or  J    (-a). 

Then  by  Parseval's  theorem,  the  right-hand- side  of 

o 

equation  (98)  becomes  zero  because  of  the  orthogonality  of 
the  E    (a)  and  J    (a) . 

Li  Li 

Therefore  equation  (98)  takes  the  final  form 

+  0O 

J       G4(a,6)c7*(a)Jz(a)da  =  0  (104) 


or 


+  0O 

I        G4(a,X/X')  |«7  (a)  |  2da  =  0  (105) 


where 

3  =  2tt/A'  . 
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So,  one  can  see  the  integrand  is  dependent  on  the  ratio 

of  A/A1  or  effective  dielectric  constant,  £     ,  defined  by 

reff 

£r     s  -C  XT  )2- 

reff     A 

Using  the  preceding  current  distributions . for  single  and 
coupled  lines,  the  phase  velocity  characteristics  can  be 
determined  by  finding  those  values  of  the  propagation  con- 
stant that  give  a  zero  value  for  the  integral  at  a  particular 
frequency. 

It  is  clear  that  once  equation  (105)  is  integrated  the 
dependence  on  the  variable  a  disappears,  so  one  can  state 
that  the  whole  process  is  really  a  function  of  frequency 
and  the  structure's  physical  characteristics,  namely,  width 
of  the  strips  W,  and  in  the  case  of  coupled  strips,  separa- 
tion S  between  the  strips,  thickness  of  substrate  D,  and 
the  dielectric's  relative  permittivity  z    . 

By  observing  the  behavior  of  each  term  equation  (105) 
one  can  conclude  integration  of  the  \J    (a) | 2  term  always 
gives  a  positive  area.   In  the  three  configuration  of 
microstrip  which  were  studied,  the  Fourier  transforms  of 
the  current  distributions  always  provided  a  multiplicative 
factor  of  (Sina/a)2  which  gives  rise  to  a  curve  having  a 
large  main  lobe  and  low  side  lobes  in  the  a-domain. 

Several  plots  of  G.(cx,A/A')  and  the  G.  (a  ,  A/A  ' )  •  |  J    (a)  |  2 
term  between  limits  of  a;  -3000<a<3000  for  different  A/A' 
values  is  shown  in  Figure  4  and  Figure  5  respectively. 
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Figure    4.      G.(a,X/X')    versus    a 
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-003*** 


■  j.       I  T  -T,    11JW— Iff 

052  x  ioJ 


Figure    5.       G4  (a , A/A ' )  | J    (a)  | 2    versus    a 
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The  family  of  curves  is  symmetrical  to  the  y-axis  and 
for  decreasing  values  of  A/A',  moves  along  the  y-axis  with- 
out losing  the  symmetry.   The  total  area  under  it  monotoni- 
cally  decreases. 

The  value  of  A/A'  which  gives  a  zero  total  area  under 

the  curve  is  the  root  of  equation  (105)  .   This  gives  the 

effective  dielectric  constant  e     and  the  phase  velocity. 

reff 

D.   THE  CHARACTERISTIC  IMPEDANCE  IN  TERMS  OF  THE  DISPERSION 
CHARACTERISTICS 

In  general,  one  can  express  the  average  power  in  terms 

of  the  strip  current  as, 

PAVP  =  \   I2  Z  (106) 

AVE    2   o   o  K        J 

z 

from  which 

2P 
AVF 
Z0  =  -yF^  ,  (107) 

o 

z 

provided  I   is  a  unique  current. 
r         o         H 

In  this  study,  it  was  assumed  in  the  dispersion  charac- 
teristic part  that  the  surface  current  across  the  strips 
were  uniform;  furthermore,  the  current  can  be  expressed,  as 


I    =  0   H-d£  =        (H    -  H   )dx  (108) 

°z  J  STRIPJ  xl  x2 


I    =/       |J  (x)|dx.  (109) 

z  JSTRIP    L 


In  the  case  of  a  single  strip, 


33 


/  I 


J.  00  =  { 


—TjT     on  the   str:  — 


0    elsewhere 


For  simplicity  I   is  arbitrarily  set  equal  to  unity. 


Thus  , 


I2   =  1 
o 

z 


(110) 


As  stated  by  Collin  [Ref .  1] ,  a  general  expression  for 
the  time  average  power  flow  in  terms  of  the  electric  and 
magnetic  fields  is 


PAVE0,y)  =  Re 


h  J  J    E  x  ttfa  da 
s 


(111) 


where 


F  x  H*. a   =  E  H*  -  E  H*. 

z    x  y    y  x 


This  equation  becomes 


PAVE(x'y)  =  Re 


(E  H*  -  E  H*)dxdy 
v  x  y    y  x^    J 


(112) 


Referring  back  to  equations  (1)  through  (6) ,  it  is  found 

that  the  equations  for  the  electric  and  magnetic  fields  can 
be  in  terms  of  Hertzian  potential  functions. 

So  substituting  equations  (3)  through  (6)  in  the  expres- 


sion for  average  power, 


PAVE(x'y)  =  h   Re 


OB-f^  "  JUHl-fy  )e-)   (-j3-^y 


9(j) 


Jp 


+  jw£:iy-  )e  J'"dxdy 


(113) 


•0  6 


H   ,JSz 


O^T7  +  W^tr  )e' 
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c-J»H 


juei|_  )e-JBzdxdy 


(113) 


Applying  Parseval's  theorem 


AVE        4tt 


Re/ 


+00      +00 


-00  -00 


n^e,         -,                 3$    (a,y) 
a3$    (a,y)    -    jooy ^    >* J 


.09wYa,y)  ,e    r        . 

33 — ^  • } <-      -    aa)e$       (a,y) 


j3-gy  7  +ouoy$    (a,y) 


it  e  1 

8a$       (a,y)    -    jcoe ^ — *-^M       dady). 


(114) 


where    the    Fourier   transform   of   34>/3x   has   been   applied    from 
equation    (114).      Performing   the    operations    indicated, 

( 

+00      +00    , 

h 


PAVE         4tt    Re 


-  00      -  00 


a2  3we  I  <f(a  ,y) 


2  „    1 3$    (a,y)     1  2 


CO 


n  1  3$  (a  ,y)  12     2   ouh,    -,  1  2 
3c| — ^  tzj-jz  -  orcoy3.|$  (a,y)|z 


(115) 


/         *e                 h 
•  i  2  Lfh ,    .3$   (a,y)    *e*,    .3$  (a,y) 
jak2  4  (a,y) ^  >/J  +  $   (a,y) ^y  i  l 


-    ja3 


h* 


eCa,7)3$,(a^  *  *h%,y)^V^) 


3y 


3y 


dady 


where 


k2  =  w2ye. 
After  obtaining  the  general  expression  for  average  power, 
one  can  apply  it  to  both  regions  for  further  specific 
equations . 


35 


In  Region  1, 


*lO,y)  =  Ae(a)e 


*J(a,y)  =  An(a)e  L 


(116) 
(117) 


3*l(a,y) 

ay 


-y^y-D) 

-  YjA  (a)e 


(118) 


3*1 (a, y)     *   h    -YiCy-D) 


(U9) 


Substituting  the  preceding  equations  in  (115)  one  can 
obtain 


+  00     +0O 


P-     =  -^  Re 
XAVE    47T 


2Yl(y-D) 


D 


-a  pe.  I A  (a)  I 


(120) 


■a)3y1Yj|Ah(a)  |  2  -  coBe-jjJlA^a)  |  2  -  a203U13  |  Ah  (a) 

dadyl 


jak2(2y1Ah(a)Ae*(a))  +  j  a32  (2Y-,Ae  (a)Ah*  (a)) 


It  is  clear  from  the  above  expression  that  one  can  inte 
grate  with  respect  to  y  to  obtain 


+  0O 


AVE 


_1 

8tt 


a2  +  y5«  ' 


Y: 


3a)£1  I  A  (a) 


(121) 

da. 


+  Ba)y1|Ah(a)  |  2  +  2a(32+k2)Re  f-  j  Ae  (a)  A*h  (a) 


In  Region  2,  clearly  there  are  two  expressions  for  aver- 
age power.   When  y?    is  real,  equations  (47)  and  (49)  apply 
as  stated  earlier. 
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CT  H(a)    and   BT   u(a)    are    zero    so, 


*|(a,y)    ■    B^(a)SinhY2y 


*£(<*, y)    =    c|j(a)CoshY2y 


(122) 
(123) 


a*?(a,y)  h 

— Ty =    Y2CH(a)SinhY2y 


(124) 


9$2(a,y) 

ay 


=    Y2BH(a)Coshy2y 


(125) 


and 


D 


'AVE 


II 


\  Re  /  ff         Sinh2Y2y[a26oae2|B^(a) 


HYP"   "0 
REGION 


+    eooy2Y^|c|j(a)|2    +  ja62Y2Be  (a)  c|}*  (a) 

-  jaY2k2CH(a)B^  (a)]dady  +  Cosh2Y2y        (lib} 

[a26wy2|c[j(oO  |2    +    Bu>e2y|  |  B®(o)  |  2 


+    jaB2Y2B^(a)C^(a)         j  ak2  Y2C*(a)  B®*  (a)  ]  dadyj) 

The  y-dependence  of  the  average  power  in  Region  2  dis 

appears  through  corresponding  integrations  as  indicated 

below: 

D 

f  Sinh2Y9D  -  2y9D 

/   Sinh2Y2ydy  = ,i-  (12 


0' 


2y. 
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D 


/   Cosh2Y2ydy  = 


Sinh2y2D  +  2y2D 
2Vo 


0 

Thus, 

^AVE 


16tt 


Re 


(Sinh2y2D 


2Y2D) 


H 


HYP 
REGION 


a  3u)£ 
Yo 


2  ,„e 
H 


:a)|2  +  3ooy,YJC?;(a)  |2 


2  '2'  H 


+  ja32B^(a)C*h(a)  -  jak2C^(a)B*e(a) 


(Sinh2Y2D  +  2y2D) 

HYP 
REGION 


(128) 


da 


a  3wy 

y7~ 


-  IcjjCcO 


+  ea)e2Y2|B{j(a)  |2  +  ja32B*  (a)  CH"(a) 


jak2CH(a)B^  (a) 


da' 


(129) 


It  is  also  clear  that  for  imaginary  y~    the  average  power 
expression  can  easily  be  obtained  with  substitution  of  (48) 
and  (50)  and  the  corresponding  trigonometric  case  electric 
and  magnetic  field  constants  i.e.  equation  (72)  and  (74). 
A  detailed  exposure  to  such  analysis   is  given  in  Appendix 
C. 

Then  for  imaginary  y?    Region  2 : 
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'AVE, 


16tt 


Re 


(2Y'2'D 


Sin2Y'2'D)  [ 


a  3we. 


'TRIG 
V.  REGION 


-ii  I  rn 


,*h 


|B^(a)|2  +  3ojy2Y2'|C^(a)  |2  +  a32B^  (a)  CT"  (a) 


[ 


ak2C!J(a)BTe(a)]da  +  /       (2y2"D  +  Sin2Y2D) 

TRIG 
REGION 

a23wu9  ,  *>, 

„   |c"(a)  |2  +  3a)e2Y,2,|B®(a)  |2  -  a32B^  (a)CTn  (a) 


,h 
"2WT 


ak2Ci!(a)B^  (a)]  da  >  . 


(130) 


Total    power    in   Region    2   becomes: 


=    P. 


+   P. 


'AVE 


'AVE, 


'AVE, 


(131) 
JH  ~V"T 

As  strip  width  decreases  one  observes  a  widening  a-domain 
power  distribution  in  Region  1.   Also  for  higher  values  of 
dielectric  constant  where  the  energy  is  more  concentrated 
into  the  substrate,  a  wider  a-domain  power  distribution  is 
obtained  in  Region  2. 

Finally,  characteristic  impedance  can  be  found  as 


2(P 


+  P. 


AVE 


'AVE 


) 


'0 


I2 
o 


(132) 


In  all  cases,  geometry  is  taken  into  account  in  the 
current  distribution  expressions. 

Plots  of  a-domain  power  distribution  for  Regions  1  and 
2  are  shown  in  Figure  6  and  Figure  7  respectively.   By 
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Figure  6.   Region  1  Average  Power  Distribution  in  a-Domain 
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Figure  7.   Region  2  Average  Power  Distribution  in  a-Domain 
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taking  the  inverse  Fourier  transforms  of  P1     and  P?     one 

iAVE      ZAVE 
could  obtain  the  x-domain  power  distribution. 

It  is  important  to  observe  and  prove  the  continuation  of 

P-     and  P?    curves  for  the  critical  points  where  y? 
■•"AVE      ZAVE  L 

switches  from  real  to  imaginary  character  or  vice  versa. 

The  enlarged  views  of  P-     and  V  given  in  Figures  8 

"""AVE      ZAVE 
and  9  clearly  show  the  smoothness  of  transition  from  one 

region  to  another. 


42 


K-5C.RLF    5.  OOF.-OI   UNITS    INCH 
y-5CRLF>l.  00: +Q1   UNJT5    ]NCH 


Figure  8.   P-     versus  aD  at  1  GHz 
"•"AVE 
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X  $CRLF-=5.00E-01  JNIT3  INCH 
Y-3C,RLF.=i.0OE+O2  JNIT3  INCH 


Figure  9.   P?    versus  aD  at  1  GHz 
ZAVE 
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III.   PERTURBATION  ANALYSIS  OF  MICROSTRIP  ON  FERRITE  SUBSTRATE 


Considered  here  is  an  extension  of  the  spectral  domain 
technique  to  the  analysis  of  microstrip  on  a  ferrite  substrate 
A  familiar  perturbation  expression  for  the  propagation  con- 
stant of  a  hollow,  closed  boundary  waveguide  with  a  change  in 
material  is  given  by  Helszajn  [Ref.  5]  as: 


(f  +  r*)  = 


jui   /jf(eo[AXe]-F-E*  +  yQ  [AXJ  -H'-H*)  da 


(E*  x  H'  +  ExH*)  «a  da 


(133) 


which  primes  denote  perturbed  quantities.   This  expression 
also  is  valid  for  open  boundary  structures  which  support 
bound  waves  [Ref.  7]  as  is  the  case  here. 

Consider  the  case  of  a  microstrip  on  ferrite  with  mag- 
netization perpendicular  to  the  substrate.   In  this  case, 
assuming  the  lossless  dielectric  substrate  problem  is  solved, 


there  is  no  change  in  electric  susceptibility  so  Ax 
The  change  in  magnetic  susceptibility  is  given  by 


=  0. 


[Ax 
Am_ 


X     Ox 
Axx     Axz 


0    0    0 


Xzx     xzz 


(134) 


where  magnetic  loss  is  permitted.   The  denominator  of  (133) 


is  just  4P.yp  so  one  obtains 


[Ax  ] -H' -H*da 


r1  +  r*  =  jooy 


11    12 

O  0 


(135) 
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Evaluation  of  the  numerator  of  (135)  leads  to: 


I"  +  r*  =  jooy 


/o 


(X   H  H*  +  X   H  H*)dydx 

^AXX  X  X    Azz  z  ZJ     7 


11   I2 

o  o 


(136) 


where  unperturbed  fields  have  been  used  as  a  first  approxi- 
mation to  perturbed  fields  and  where  because  of  the  geometry 
no  demagnetization  of  r.f.  fields  takes  place.   No  cross 
terms  involving  products  of  H   and  H   appear  in  the  numera- 

Lt  A. 

tor  of  (136)  because  this  product  is  odd  and 


+  00 


+  oo 


H  H*dx 
z  x 


H  H*dx  =  0 
X  z 


(137) 


Referring  back  to  equations  (2)  and  (5) ,  substituting 
in  equation  (156)  and  applying  the  Parseval's  theorem  one 
obtains: 


r '  +  r*  = 


3toy 


+  oo   D 


2Z  I2    2  7T   /  } 

°  °Z         -«    0 


'XX 


±j3(-jo*2(a,y) 


3$2Ca,y) 

+  ja)£0 k 

J       2         3y 


3*|(a,y) 


+  j3(-ja$2(a,y)  -  jwe2 g- ) 


+  xzzkc  ^2^a,y^$2  (a>y)}    d>rda 


(138) 


where 


Azz   Axx 


This  propagation  constant  expression  holds  for  ±  z 
directed  traveling  waves  depending  upon  the  choice  of  upper 
or  lower  sign  in  (138). 
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One  might  attempt  to  solve  this  equation  by  defining 
regions  of  a  where  y?  has  imaginary  and  real  values  which 
are  indicated  by  equations  (36)  and  (37) . 

Also  substituting  previously  calculated  corresponding 

h  e 

equations  of  ^-(ajy)  and  (8<I>  (a,y)  )/9y  for  hyperbolic  and 

trigonometric  cases  into  (138)  and  integrating  out  the  y 

dependence  the  following  is  obtained: 


r«   +  r 


JO)U 


16tt 


Z    I2 
o   o 


z     -/TRIG 
REGION 


2y2D    +    Sin2y2D 

7:      / 


Xxxf("B2a2    +   k*    )|c!Jl(a)|2    t   joj^Bcr^' 


hi  2  I  De 


(-jC;i(a)B,r"(a)    -    jC^a)B^(a))    -    (we^Y")2  I  BX  (a) 

11  1  x  £m       it  1 


da    + 


jam 


Sinh2Y?D 


2Y2D 


16tt    Z    I2 
o   o 


Y- 


HYP 
REGION 


i 


Xxx[("62a2+kJ  )  |c|)(a)  |2±jaJ£26aY2(C^(a)B{ei:fe(a) 


c|j*(a)B®(a))  -  (a)£:2Y2)2|B^(a)|2]\da.    (139) 


Evaluating  the  difference  between  perturbed  and  inper 
turbed  propagation  constants, 


r'  +  r*  =  Aa  +  j A3 

r«  +  r*  =  a'  +  j3'  -  j3 


(140) 


where 
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Act  =  a'  and  A3  =  3'  -  3 


Also,  as  stated  in  Reference  5 


xxx   xxx    Jxxx 


(141) 


Substituting  (140)  and  (141)  into  equation  (139)  and 
separating  real  and  imaginary  parts  gives  expressions  for 
a'  and  3',  loss  and  phase  shift  with  ferrite  respectively 


Define  a  common  expression 


K 


jam 


16tt  '  Z  IM 
o  o 


2Y2  +  Sin2Y'2'D 


TRIG 
\,  REGION 


v  rr 


2'  T 


-  (a)C2Y^)2|B^(a) 


da  + 


'T 


Sinh2y2D 


2y2D 


HYP 
REGION 


(-32a2+k^  )|c[](a)|2  ±  2coe23aY2(-jC*h(a)B^(a)) 
2  \ 


(we2Y2)2 lBn(a) 


da 


Then 


Aa+  =  K+X^x 
AP+  =  K+X'x 
where  x'   and  x"   are  defined  in  Reference  5. 

AA  A.  A. 


(142) 

(143) 
(144) 


48 


Loss  in  db  per  wavelength  versus  frequency  and  $.p/3 
versus  frequency  are  plotted  in  Figures  10  and  11. 
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Figure  10.   Loss  per  wavelength  in  db  versus  frequency  (GHz) 


50 


8 


s 


> 


oo: 


QC5 


010 


us 


023 


020 


K-SCflU>5.00E-+00  units  inch. 
V-SCRLE----2.00E-01   UNITS    INCH, 

Figure  11.   6^/3  versus  frequency  (GHz) 
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IV.   COMPUTER  PROGRAMMING 

A  computer  program  is  written  in  FORTRAN  IV  language 
that  first  calculates  the  phase  velocity  or  A/A'  ratio  for 
a  given  geometry  of  microstrip  and  parameters  of  substrate. 
Then  the  program  proceeds  and  calculates  the  characteristic 
impedance  for  the  same  set  of  data  using  results  found  in 
the  first  phase  of  the  calculation  namely  A/A'.   If  the  sub- 
strate is  ferrite,  also  phase  shift  and  loss  per  wavelength 
are  calculated  and  curves  are  generated. 

The  program  consists  of  one  main  routine  and  three  sub- 
routines.  The  main  routine  starts  by  calculating  all  con- 
stants and  parameters  which  arc  needed  tc  evaluate  the 
various  integrals.   First  equation  (105)  is  integrated  and 
evaluated  for  an  arbitrary  initial  A/A'  value.   This  process 
continues  for  different  values  of  A/A'  until  zero  is  ob- 
tained as  the  value  of  integral.   This  computer  time  consum- 
ing iteration  process  is  shortened  by  a  simple  root  finding 
routine  that  determines  the  correct  region  of  A/A'  and  elim- 
inates others.   On  the  average,  after  9  iterations  equation 
(105)  gives  a  value  which  will  yield  zero  value  for  the  in- 
tegral with  an  accuracy  of  0.011.   The  exact  number  of  steps 
will  change  depending  on  desired  accuracy  and  physical 
parameters . 

Evaluation  of  the  integral  is  done  by  standard  IBM 
360/67  system  library  subroutine  DQG24  which  uses  a  24 
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point  Gaussian  Quadrature  formula  which  integrates  poly- 
nomials up  to  47th  degree. 

This  routine  has  been  found  to  accurately  integrate  the 
type  of  curves  encountered  in  this  study. 

The  DQG24  routine  requires  an  external  function  sub- 
program which  defines  the  polynomial  to  be  integrated  for 
imaginary  and  real  values  of  y~.      Two  different  subprograms 
GZR,  GZIM  were  used  for  this  purpose  and  to  test  and  set 
the  limits  of  integration,  for  real  and  imaginary  character 
of  Yo  respectively. 

It  is  also  proven  that  increasing  the  limits  of  inte- 
gration or  using  a  32  point  integration  routine  which 
evaluates  polynomials  up  to  64th  degree  does  not  change 
the  result  more  than  0.1% . 

Another  subroutine  called  TRAN  supplies  the  expression 
for  the  transformed  current  distribution  depending  upon  the 
physical  configuration  of  the  microstrip. 

As  far  as  the  theoretical  development  and  programming 
is  concerned  the  only  approximations  used  are  the  repre- 
sentation of  the  axial  surface  current  density  by  a  square 
pulse  and  the  assumption  that  the  transverse  surface  cur- 
rent density  is  zero.   Different  forms  of  current  distribu- 
tions were  tried  such  as 

f 
coshx   on  the  strip 


J  (x)  = 


elsewhere . 
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It  was  concluded  that  due  to  the  many  zero  crossings  of 
this  transformed  current  distribution  curve  and  the  greater 
number  of  side  lobes  which  adversely  affected  programming 
efficiency  its  use  was  not  justified  even  though  it  more 
accurately  approximates  the  true  current  distribution. 

Characteristic  impedance  was  calculated  as  a  next  step 
by  evaluating  equations  (129)  ,  (130)  and  (132)  using  dif- 
ferent parts  of  the  same  routines  used  in  the  first  part 
of  the  calculation. 

If  the  substrate  is  ferrite  the  program  proceeds  to 
another  section  which  evaluates  equations  (142)  through 
(144)  and  calculates  loss  and  phase  shift  values  as  a  func- 
tion of  operation  frequency. 

It  may  be  pertinent  to  quote  the  typical  computational 
time  using  the  method  outlined  in  this  study. 

From  the  GO  step  the  average  time  for  calculating  X/X 


Z  ,  and  loss  and  phase  shift  (if  necessary)  for  a  single 
frequency  was  between  1.5  to  1.9  sec  on  an  IBM  360/67. 
The  computer  program  accepts  the  following  data: 

-  Substrate  parameters,  e   ,  u 

r2    r2 

-  Thickness  of  substrate,  D,  in  millimeters. 

-  The  separation  between  conductors,  S,  in  millimeters 
(only  for  coupled  strip) . 

-  Initial  frequency  (GHz) . 

The  program  starts  calculating  values  for  the  initial 
frequency  and  repeats  the  procedure  by  increasing  the  fre- 
quency in  1  GHz  steps  up  to  50  times.   Means  are  also  pro- 
vided to  compare  the  half  wavelength,  X/2,  with  the  thickness 
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of  the  substrate,  D.   For 


D  = 


the  program  stops  due  to  the  presence  of  higher  order  modes 
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V.   NUMERICAL  RESULTS  AND  COMPARISONS 

The  dispersion,  characteristic  of  the  microstrip  as  a 
function  of  frequency  has  been  evaluated  for  a  variety  of 
geometries  and  these  results  are  presented  both  in  tabular 
and  plotted  form.   In  addition,  characteristic  impedance  as 
a  function  of  frequency  is  presented  for  the  same  geometries 
and  dielectric  constant. 

Figures  12  through  15  show  the  frequency  variation  of 
wavelength  and  characteristic  impedance  of  a  single  strip 
for  various  values  of  W/D.   Similar  families  of  curves  for 
coupled  strips  appear  in  Figures  16  through  19.   In  this 
case  the  parameter  is  S/D  and  W/D  is  held  constant. 

The  results  indicate  that  higher  frequencies  and  wider 

strips  increase  the  effective  dielectric  constant,  e     , 

reff 
and  the  characteristic  impedance  Z   of  both  single  and 

coupled  strips.   The  spacing  S  between  coupled  strips  also 

has  an  effect  on  e     .   Increasing  spacing,  S,  increases 

reff 
e     for  odd  and  even  modes.   The  effect  on  characteristic 

reff 
impedance  depends  upon  the  mode  of  propagation.   Increasing 


spacing,  S,  increases  Z   for  the  odd  coupled  mode  but  causes 

a  decrease  in  Z   for  the  even  coupled  mode. 

o  r 

All  available  published  data  is  based  on  quasi-static 
analyses  except  in  a  few  cases.   Basically,  in  this  study 
a  third  dimension,  frequency  or  D/X,  is  introduced.   Two 
3-D  plots  are  given  in  Figures  20  and  21  to  illustrate  this 
point.   The  x-axis  is  D/X,  the  y  axis  W/D,  and  the  z-axis 
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Figure  12.   Single  Strip  X'/X  versus  D/X  curves  e    =  6.0 
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Figure  13.   Single  strip  Z   versus  D/X  curves  e. 
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Figure  14.   Single  strip  X'/X  versus  D/A  curves  e 
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Figure  15.   Single  strip  Z   versus  D/A  curves  e    =  12.0 
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Figure    16.      Coupled    strip    even    and   odd   modes    A'A    versus    D/X 
curves    e        =    6 . 0   W/D   =    1 . 54 . 
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Figure  17.   Coupled  strip  even  and  odd  modes  Z   versus  D/A 
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Figure  18. 


Coupled  strip  even  and  odd  modes  A '/A  versus 
D/A  curves  e    =12.0  ltf/D  =  0.8 
r2 
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Figure  19.   Coupled  strip  even  and  odd  modes  Z   versus  D/.\ 
curves  e    =12.0  W/D  =0.8.       ° 
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Figure  20.   Three  dimensional  plot  A'A,  W/D,  D/X  single  strip 
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Figure  21.   Three  dimensional  plot  Z  ,  W/D,  D/X  single  strip 
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X'/X  and  Z   respectively.   One  can  find  the  X' /X  versus  W/D 
and  Z   versus  W/D  surfaces  already  plotted  in  different 
references,  i.e.  [Ref.  10]  but  only  for  zero  frequency. 

This  newly  added  dimension  makes  it  possible  to  find 
data  for  a  specific  frequency  in  addition  to  a  certain 
geometry. 

In  Table  I  the  values  found  for  characteristic  impedance 
with  two  different  methods  in  Reference  3  and  generated  from 
this  study  are  compared.   This  comparison  reflects  the  agree- 
ment in  the  low  frequency  range. 

Also  some  of  the  tabulated  data  in  Reference  11  is  com- 
pared with  data  obtained  from  this  study  in  Table  II  for 
coupled  strips. 

One  of  the  few  studies  of  frequency  dependence  was  by 
Haddad  [Ref.  9].   Table  III  shows  data  taken  from  Haddad's 
work  and  from  this  study  which  clearly  indicates  the  agree- 
ment also  for  low  frequency  ranges.   Although  Haddad's 
analysis  was  frequency  dependent,  the  data  presented  by  him 
was  restricted  to  relatively  low  frequencies. 

Finally  three  sample  computer  output  pages  are  given 
in  Tables  IV,  V  and  VI . 

Although  the  wavelength  and  characteristic  impedance  of 
microstrips  on  dielectric  substrate  depend  only  upon  D/X, 
the  normalized  frequency,  this  is  not  true  for  ferrite  sub- 
strates.  Since  the  computer  program  will  calculate  results 
for  either  case  it  is  necessary  to  provide  actual  dimensions, 
D,  W,  S  etc.  as  inputs.   Where  both  frequency  and  D/X  are 
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TABLE  I.   Characteristic  Impedances  of  Microstrip  Transmission 
Lines  Single  Strip. 


%'- 

6.0 

W/D 

Z*(ohm) 

Z+Cohm) 
o 

Z^(ohm) 

0.1 

135.455 

134.352 

137.64 

0.2 

113.272 

112.255 

115.41 

0.4 

91.172 

89.909 

93.25 

0.7 

73.613 

71.995 

73.53 

1.0 

62.713 

60.970 

64.45 

2.0 

43.149 

41.510 

44.15 

4.0 

27.301 

26.027 

27.23 

0.0 

13.341 

12.485 

12.58 

r2 


=  16.0 


0.1        85.9659  87.762  87.70 

0.2         71.6954  73.025  73.43 

0.4  '       57.4999  58.110  59.20 

0.7        46.2344  46.217  47.83 

1.0         39.2512  38.948  40.72 

2.0         26.7555  26.248  27.73 

4.0        16.7210  16.300  16.96 

10.0  8.0385  7.807  7.77 

*  Characteristic  impedance  obtained  by  method  of  moment  [Ref.  3] 

f  Characteristic  impedance  obtained  by  conformal  mapping  [Ref.  3] 

s  Characteristic  impedance  obtained  by  spectral  domain 
method  in  this  study. 
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Table  III.   Frequency  Dependence  of  Single  and  Coupled 
Microstrip . 


SINGLE 
W/D  =1.0 

e*       et 

reff     ref£ 

1GHZ   6.78      6.78 
2GHZ   7.05      6.98 


z* 

0 

zt 

0 

1GHZ  48.5ft 

50.88ft 

2GHZ  50.0ft 

51.22ft 

COUPLED  COUPLED 

(EVEN  MODE)  (ODD  MODE) 

W/D  =1.0  S/D  =0.4   W/D  =1.0  S/D  =0.4 


xeff 
7.39 
7.8 


Z* 
o 


61.8ft 
63.0ft 


et 

reff 

7.43 
7.72 

Zt 
o 

62.24ft 

62.8  ft 


eff 

5.84 
5.87 


Z* 
o 


;t 

re££ 
5.85 
5.94 


Zt 
o 


36.2ft        39.31ft 
36.  2ft        39 . 18ft 


*  Haddad  results 

t  Spectral  Domain  Transform  method  used  in  this  study 
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Table  IV.   Sample  Computer  Output  for  Single  Strip 


W  =  1mm 


D  =  2mm 


e    =12    W/D  =  0.5 
r2 


FREQUENCY 


D/LAMBDA 


L/LPR 


LPR/L 


CH  IMPEDANCE 


1.000  GHZ 

2.000  GHZ 

3.000  GHZ 

4.000  GHZ 

5.000  GHZ 

6.000  GHZ 

7.000  GHZ 

3.000  GHZ 

9.000  GHZ 

10.000  GHZ 

11.000  GHZ 

12.000  GHZ 

13.000  GHZ 

14.000  GHZ 

15.000  GHZ 

16.000  GHZ 

17.000  GHZ 

18.000  GHZ 

19.000  GHZ 

20.000  GHZ 


0.00667 
0.01333 
0.02000 
0.02667 
0.03333 
0.04000 
0.04667 

0.06000 
0.06667 
0.07333 
0.08000 
0.08667 
0.09333 
0.10000 
0.10667 
0.11333 
0.12000 
0.12667 
0.13333 


2.75971 
2.78140 
2.80726 
2.83508 
2.86373 
2.89254 
2.92108 
2.94907 
2.97628 
3.00258 
3.02787 
3.05206 
3.07512 
3.09702 
3.11775 
3.13733 
3.15577 
3.17311 
3.18939 
3.20464 


0.36236    62.59  OHM 


0.35953    62.70  OHM 


0.35622    63.07  OHM 


0.35272    63.72  OHM 


0.34919    64.68  OHM 


0.34572    65.96  OHM 


0.34234    67.56  OHM 


0.33909    69.51  OHM 


0.33599    71.81  OHM 


0.33305    74.47  OHM 


0.33027    77.51  OHM 


0.32765    80.94  OHM 


0.32519    84.77  OHM 


0.32289    89.01  OHM 


0.32074    93.66  OHM 


0.31874    98.73  OHM 


0.31688   1C4.24  OHM 


0.31515   110.17  OHM 


0.31354   116.52  OHM 


0.31205   123.31  OHM 
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Table  V.   Sample  Computer  Output  for  Coupled  Strip  (Odd  Mode) 

W  =  1.6mm   D  =  2.0mm   S  =  0.8mm   W/D  =0.8    S/D  =0.4   e   =12 

r2 


FREQUENCY 


D/LAMBDA 


L/LPR 


LPR/L 


CH  IMPEDANCE 


1.000  GHZ 

2.000  GHZ 

3.000  GHZ 

4.000  GHZ 

5.000  GHZ 

6.000  GHZ 

7.000  GHZ 

8.000  GHZ 

9.000  GHZ 

10.000  GHZ 

11.000  GHZ 

12.000  GHZ 

13.000  GHZ 

14.000  GHZ 

15.000  GHZ 

16.000  GHZ 

17.000  GHZ 

18.000  GHZ 

19.000  GHZ 

20.000  GHZ 


0.00667 
0.01333 
0.02000 
0.02667 
0.03333 
0.04000 
0.04667 
0.05333 
0.06000 
0.06667 
0.07333 
0.08000 
0.08667 
0.09333 
0.10000 
0.10667 
0.11333 
0.12000 
0.12667 
0.13333 


2.61034 
2.61858 
2.63141 
2.64815 
2.66821 
2.69112 
2.71640 
2.74366 
2.77249 
2.80248 
2.83325 
2.86441 
2.89561 
2.92649 
2.95675 
2.98613 
3.01442 
3.04146 
3.06714 
3.09139 


0.38309    39.05  OHM 


0.38189    38.97  OHM 


0.38002    38.89  OHM 


0.37762    38.84  OHM 


0.37478    38.87  OHM 


0.37159    39.00  OHM 


0.36813    39.26  OHM 


0.3644c    39.69  OHM 


0.36069    40.31  OHM 


0.35683    41.16  OHM 


0.35295    42.26  OHM 


0.34911    43.65  OHM 


0.34535    45.34  OHM 


0.34171    47.37  OHM 


0.33821    49.75  OHM 


0.33488    52.51  OHM 


0.33174    55.67  OHM 


0.32879    59.23  OHM 


0.32604    63.22  OHM 


0.32348    67.65  OHM 
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Table  VI.   Sample  Computer  Output  for  Coupled  Strip  (Even  Mode) 


UNI. 6mm  D=2.0mm  S=0.8mm   W/D  =  0.8   S/D  =  0.4   e   =12 

r2 


FREQUENCY     D/LAMBOA 


L/LPR 


LPR/L 


CH  IMPEDANCE 


1.000  GHZ 

2.000  GHZ 

3.000  GHZ 

4.000  GHZ 

5.000  GHZ 

6.000  GHZ 

7.000  GHZ 

8.000  GHZ 

9.000  GHZ 

10.000  GHZ 

11.000  GHZ 

12.000  GHZ 

13.000  GHZ 

14.000  GHZ 

15.000  GHZ 

16.000  GHZ 

17.000  GHZ 

18.000  GHZ 

19.000  GHZ 

20.000  GHZ 


0.00667 
0.01333 
0.02000 
0.02667 
0.03333 
0.04000 
0.04667 
0.05333 
0.06000 
0.06667 
0.07333 
0.08000 
0.08667 
0.09333 
0.10000 
0.10667 
0.11333 
C. 12000 
0.12667 
0.13333 


2.92702 
2.96430 
3.00492 
3.04470 
3.08186 
3.11570 
3.14609 
3.17316 
3.19719 
3.21849 
3.23740 
3.25422 
3.26921 
3.28263 
3.29467 
3.30552 
3.31533 
3.32423 
3.33234 
3.33975 


0.34164    64.22  OHM 


0.33735    64.38  OHM 


0.33279    64.92  OHM 


0.32844    65.82  OHM 


0.32446    67.00  OHM 


0.32096    68.42  OHM 


0.31786    70.00  OHM 


0.31514    71.70  OHM 


0.31276    73.49  OHM 


0.31070    75.35  OHM 


0.30889    77.26  OHM 


0.30729    79.20  OHM 


0.30588    81.16  OHM 


0.30463    83.14  OHM 


0.30352    65.13  OHM 


0.30252    87.13  OHM 


0.30163    89.14  OHM 


0.30082    91.15  OHM 


0.30009    93.17  OHM 


0.29942    95.18  OHM 
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shown  as  printed  output  the  calculation  is  based  on  the  sub 
strate  thickness,  D,  provided  as  an  input.   A  width  of  D  = 
2mm  has  been  used  in  the  calculations  presented  here. 
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VI.   CONCLUSIONS 

In  this  thesis  a  hybrid  mode  analysis  of  single  and 
coupled  microstrips  (of  equal  width)  on  dielectric  substrate 
was  described.   The  resulting  equations  were  solved  by  the 
method  of  moments  applied  in  the  spectral  or  (Fourier)  trans- 
form domain.   A  computer  program  was  developed  to  calculate 
wavelength  and  characteristic  impedance  versus  frequency. 
The  original  results  obtained  by  this  method  were  compared 
in  the  low  frequency  limit  with  the  quasi-static  results 
of  several  other  investigators  and  were  found  to  be  in  ex- 
cellent agreement. 

The  method  was  extended  to  microstrip  on  ferrite  where 
the  propagation  constant  was  calculated  using  perturbation 
theory.   Numerical  evaluation  of  resulting  expressions  is 
again  carried  out  in  the  spectral  domain.   A  computer  program 
has  been  developed  for  this  analysis  also. 
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APPENDIX  A 
AUXILIARY  VECTOR  POTENTIAL  FUNCTIONS 

A.   MAGNETIC  HERTZIAN  VECTOR  POTENTIAL  FUNCTION 

Consider  a  homogeneous,  source-free,  isotropic  region; 
hence,  there  is  no  charge  density  and 

V«E  =  0  (Al) 

From  Vector  Algebra,  it  is  known  that  the  divergence  of 
a  vector  is  zero  if  the  vector  is,  in  turn,  the  curl  of 
another  vector;  therefore,  one  can  state, 

V«Vx?h  =  0    >f  ?h  (A2) 

E  =  -jaj(aVX7Tr  .  (A3) 

For  time-varying  fields 

VxH  =  jue  E  (A4) 

so  applying  equation  (A3) , 

VxH     =     J03£  (-  JOOUVXTT,    ) 

=    k2Vx7r, 


=   k2Vx(7T,     +   k"    V<j>) 


(A5) 
Vx(tt,     +   k    ~V<J>) 

H   =    k2?,     +   V<J>     .  (A6) 


where , 

k2     =    0)2U£. 

Similarly,    for   time-varying   fields, 

VxE    =    -jcoyiT  (A7) 
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so,  applying  equation  (A3), 

Vx(-jcoyVx7h)    =    -jojy(k2¥h   +    V<J>)  (A8) 

W-7,     -    V27,     =   k2?h   +   V«f>.  (A9) 

Up  to  this  point  <$>   has  been  defined;  arbitrarily  choose 
<J>  =  V«?h.  (A10) 

Then,  equation  (A9)  becomes, 

V2Tfh  +  k2^h   =    0.  (All) 

Also,    it    is   known    that 

V»B    =    yV»H   =    0 

uV«  (k2%    +    Vcf>)    =    0  '  (A12) 


k2V»7fh    +    V'Vcf)    =    0. 

Summarizing,    for  TM  modes, 

E    =    -jojuVxtF,  (A13) 

H   =    k2?h   +    W»?h  (A14) 

-     VxVXTT,    . 

h 
Similarly,    for   TE   modes, 

H   =    jojeVx?h  (A15) 


E     =     k^TT         +     VV'TT 


=     VXVXTT      . 

e 


(A16) 


B.      TE   AND   TM  MODES    FROM  VECTOR   POTENTIALS 
1 .      TE   Modes 
E    =    -ju>yVx7h  (A17) 
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H"  =  k2Fh  +  W-7rh 
c 

=■  =•  -h  (A18) 

=     VxVXTT  V  J 

where    tt      satisfies: 

V2¥h   +   k2?h  =    0  (A19) 

c  J 

For   TE   modes,    E      =0:    therefore 
~   h        ~    h 

Ez  =  -^  t  -£  -  -57  >  fA201 

from  where, 

irjj  *   ir£  =   0.  (A21) 

Let: 

7h   =    (J)he"yza  (A22) 

Then,    from   equation    (A13) 

V2cf>h   +   k2c{)h   =    0  (A23) 

where , 

k2    =    y2    +   k2.  (A24) 

Summarizing , 

Hz    =    k2<J>h(x,y)e±YZ  (A25) 

H   =  ±ye±YZVt(})h  (A26) 

E.    =    ±i^  a"  xFT  (A27) 

t  y         Z      t  *■         ' 

2.       TM  Modes 


E    =    k2;?T      +    VV-7T  (A28) 

=     VxVXTT 

H   =    jweVx?6  (A29) 
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For   TM  Modes,    H      =0;    therefore 
^e   -   *e(x,y)e"YZIz  (A30) 


and, 


Vi<}>      +   k2cj)      =    0  (A31) 

t  c  ■ 


where , 

k2    =   y2    +   k2.  (A32) 

Summarizing, 

Ez   =    k24>ee±YZ  (A33) 

E      =    ±YVt<J>ee±YZ  (A34) 

H      =   i   2®L  J  xf  (A35) 

t                y         z      t  v         ' 

Equations  (1)  and  (2)  correspond  to  equations  (A33) 
and  (A25)  ,  respectively. 
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APPENDIX  B 
TRANSVERSE  ELECTRIC  AND  MAGNETIC  FIELDS 

Consider  Maxwell's  point  form  of  Amere's  law, 

VxH  =  jweE  (Bl) 

3H     3H  3H     9H  3H     3H 

—    c      z  y   ^        —      r        x  z    ■,        —      r        y  x. 

V^Sy  -     "SIT  3  +  ay  C  -TY  -  -5J  )  +  az   (  -gi  -  -jy   ) 

■  J^^x^  +  Vy  +  EzV  <B2> 


from  where, 


3H    3H 

z y  , 


9y    dz        J    x 
3H     3H 


=  jweEv  (B3) 


jcoeE  (B4) 


9z     9x    J    y 

9H     9H 

~£   '    -37  "  >£Ez  CB5> 

V  z 

and,  due  to  the  z-dependence  of  the  fields,  i.e,,  e'  ,  one 
can  apply 

af  e*2  =  Ye^z  (B6) 

obtaining    from   equations    (B3)    through    (B5) 

9H 

YH      =    ja)£Ev  (B7) 


9y  r»y         j—  x 

YHx    "    if   =    >£Ey  (B8) 

3H  3H 

-5^   -    -^-  =    jweE    .  (B9) 

3x3yJz  K      J 

80 


Similarly,  starting  from  Maxwell's  point  form  of  Lenz' 
law, 

VxE  =  -jtouH  (BIO) 

3E     3E  3E     3E  3E     3E 

x  k   3y     dz    J  y  3z     3x  '     z  v   3x     3y 


=  -jwyfH  a   +  H  a   +  H  a  )  (Bll) 

JV-xx    yy    z  z'  v   J 


from  where , 


3EZ 

3y 

YEy   =    -jooyHx 

3E 

^Ex  " 

=    -jooyH 

x                 y 

3Ey  . 

BE 

X                .         u 

— = —    =     -IdluH 

3x     3y 


(B12) 
(B13) 

(B14) 


Substituting  the  value  of  E   from  equation  (B13)  into 

.A. 

equation  (B7) ,  one  obtains 


3H  .      3E 

z     „  nooe  r         z  •   o  -\ 

— s—  -  yH  -1 (  — 7T—  -  jwyH  1 

3y    '  y  y   *■  3x  J   y  ** 

3H         3E 

Hy(Y2+o32ye)   =  y  ~af  -  >£  -gf  (B15) 

-     3H        3E 

H   =  XT     (Y  —it—  -  JU)£-^—  )  . 

y   k2  w   3y   J         dx    J 

Similarly,  substituting  the  values  of  E   from  equation 

(B12)  into  equation  (B8) ,  of  H   from  equation  (B8)  into 

equation  (B12)  and  of  H   from  equation  (B7)  into  equation 

(B13) ,  one  obtains 

-  For  H  , 
x' 
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3H  1         3E 

yH —  «   jwe  -   (   -5^-  +    jwyH   ) 

1    x  3x         J         Y  3y  X-* 

3H  3E 

Hx(Y2+032ye)    =   Y  _|  +   ju>e  -^  (B16) 

1  3H  3E 

x        k2     v         3x        J  3y    '. 

c  7 

-  For  E    , 

y 

3E  1         8Hz 

—^ —   -    yE      =    -iwy   —    (    — k —  +    jcoeE    ) 
3y         '    y  J         y  3x        J         y' 

3E  3H 

Ey(Y2+o)2y£)    =    Y   -gy   +   j  uy   -g|  (B17) 

-  9E  3H 

Ey  "  EJ  C     -ay  +  5«v  "if  ) . 

-  For   Ex, 

3E  3H 

Y"Ex  "  Y  "si  =   "jajy   (  "if  "  ja)£Ex) 

3E  3H 

Ex(Y2+a32y£)    -   Y  —|  -    jcoy  — i  (B18) 

1  3E  3H 

E      =    t~5-    (      — ^   -    iwy    — s—  )  . 
x        k2     v         3x         J  3y 

c  7 

Summarizing, 

1     3E         3H 

E   =  ITT  (Y  -r^  "  jwy  -s^  )  (B19) 

x    k2  k '   3x    J     3y  '  *"    ^ 

c  1 

1     3E         3H 
Ey  "  W    CY  -^  *  j»P  -jf  )  (B20) 

-,     3H         3E 

H      =    r4    (Y   -t5-  +    jwe    — 5^-   )  (B21) 

xk^v3xJ  9y 

c  7 

,  3H  3E 

H  ,    =    vr    (Y   -^   -    ju>e   -^      ).  (B22) 


y         k2     ^       3y 


3x 
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Equations  (B3)  through  (B6)  correspond  to  equations 
(B19)  through  (B22) ,  respectively. 
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APPENDIX  C 

AVERAGE  POWER  EXPRESSION  IN  REGION  2 

(TRIGONOMETRIC  CASE) 

Electric  and  Magnetic  field  potential  functions  and 
derivatives  with  respect  to  y  become 


*|(o,y)  *=  jB^Ca)SinY^'y 
*5(a,y)  =  C^(a)CosY'2'y 


(CI) 
(C2) 


— Sj- =  JY^B|(a)CosY'2,y 


(C3) 


3*5Ca.y) 


3y "  -Y2Gr(a)SinY2'y 


Substituting  into  (115) 

D 


'AVE, 


it-  Re  J 


V 


TRIG 
REGION 


II     [  [  -a23we2  I  bt  (a)SinY'-?'y  |  -wBu- 
J0    1 


(C4) 


|-Y2'C^(a)SinY2'yr-w3e2|  JY2B^(a)CosY,2y|2-a2ajy23|C^(a)CosY2*y|  2] 
•fjak2[(C^(a)CosY2'y)(-JY2B^(a)CosY2,y)    +    (- j  B^*  (a)  SinY2'y) 
(-Y^tcOSinY^y)]     ■    j  «32  [  (j  B^  (a)  SinY'2'y)  (-Y'2'C^  (a)  SinY2'y) 


,h* 


,IITj6 


+    (C^    (a)CosY2'y)(JY2B^(a)CosY^y)]]dyda 
Then, 


(C5) 
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D 


'AVE, 


4^Re<  ff         Sin2Y'2'ydy[-a2ea)e2|B^(a) 

TRTH  J  J() 


TRIG '  J0 
REGION 


"r1^  r  ~  \  De 


e     h* 


-a)3yY2,2|C^(a)|2-a2k2Y^C^(a)Bj  (a) -a32Y2'B^  (a)  C^1  (a)]  da 


D 


tt  2  I  -0^  ,-  -v  12   ^2, 


Cos2Y2ydy  [-  w6e2Y22  |B^(a)  |  2 -a2u)y2B  |  C^  (a) 


TRIG  ^0 
REGION 


ak2Y'2'c!j(a)B^*(a)  +  aB2Y2'B®  (a)c!j*  (a)  ]  da\ 


2  T 


(C6) 


Also  the  integration  with  respect  to  y  in  Region  2  may  be 
accomplished  analytically  as  indicated  below: 


D 


Sin2Y'?'ydy 


2v"D  -  Sin2v"D 

•  2 '  2 

4Y'2' 


(C7) 


.D 


Cos  Y2'ydy  = 


2y2'D  +  Sin2Y2'D 


4y'2' 


(C8) 


By  substituting  (C7)  and  (C8)  in  (C6)  and  rearranging 


one  can  easily  obtain 


Re 


'AVE 


T 


16tt    1V"   / 


(2y2'D    -    Sin2Y2fD)[ 


a    3wc 


y: 


-    |B^(a) 


'TRIG 
REGION 


3a)y9Y'9'  |c!j(a)  |  2+ak2c!j  (a)  B®*  (a)    +    a32B®  (a)  dj*  (a)  ]  da    + 


2'2'    T 


*T 


(2y2'D    +    Sin2Y2D)[ 


a    3wu 


TRIG 
REGION 


Y' 


?■      |c!f(a)  |2    +    6a)e2Y,2'|B®(a) 


2De 


aezB^(a)C1u(a)     -      k2C^(a)B^    (a)]da>. 
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(C9) 


The  final  equation  (C9)  corresponds  to  equation  (130) 
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APPENDIX  D 

C  THIS  PROGRAM  IS  DEVELOPED  TO  CALCULATE  THE  EFFECTIVE 

C  WAVELENGTH  AND  CHARACTERISTIC  IMPEDANCE  OF  SINGLE  OR 

C  CCUPLED   MICROSTRIP  ON  A  DIELECTRIC  OR  FERRITE 

C  SUBSTRATE. IF  THE  SUBSTRATE  IS  FERRITE  LOSS  PER 

C  WAVELENGTH  IN  D3  AND  NORMALIZED  PHASE  SHIFT  ARE 

C  ALSO  CALCULATED. 

C  PROGRAM  ACCEPTS  FOLLOWING  INPUT   DATA 

C  FIRST  CARD  -IDENTIFICATION   DIGITS  FORMATU3) 

C  THIS  THREE  DIGIT  INTEGER  SUPPLIES  NECCESARY 

C  INFORMATION  TO  THE  PROGRAM 

C  FIRST   DIGIT    1   MEANS  DIELECTRIC  SUBSTRATE 

C  FIRST   DIGIT    2   MEANS  FERRITE     SUBSTRATE 

C  SECOND  DIGIT    1   MEANS  SINGLE      STRIP 

C  SECOND  DIGIT    2   MEANS   COUPLED    STRIP 

C  THIRD  DIGIT     1   MEANS   EVEN   MODE 

C  THIRD  DIGIT     2   MEANS   ODD    MODE 

C  THIRD   DIGIT   IS   USED  ONLY  FOR  CCUPLED 

C  STRIP  WHEN  SECOND  DIGIT  IS  EQUAL  TO  2. ALL 

C  OTHER  TIMES  SHOULD  BE  SET  TO  ZERO. 

C  EXAMPLES 

C  110    SINGLE  STRIP  DIELECTRIC  SUBSTRATE 

C  211    COUPLED  STRIP  EVEN  MODE  ON  FERRITE  SUBSTRATE 

C  SECOND   CARD  -EPR1 , MR  1 , MR2  FORMAT ( 3F10. 5 ) 

C  EPR1  IS  THE  RELATIVE  DIELECTRIC  PERMITTIVITY  REGION  1 

C  MR1   IS  THE  RELATIVE  MAGNETIC  PERMEABILITY  REGION  1 

C  MR2   IS  THE  RELATIVE  MAGNETIC  PERMEARILITY  REGION  2 

C  THIRD  CARD-  DD,  WW1,SS    FORMAT ( 3F10 . 5 ) 

C  DD-  THICKNESS  OF  THE  SUBSTRATE   IN  MILLIMETERS 

C  W^l-WIDTH  OF  THE  STRIPS  IN  MILLIMETERS 

C  SS-  SEPERATICN  BETWEEN  STRIPS    IN  MILLIMETERS 

C  pruPTH  taro-  FREor.h   FQkMAT  (  Fi  0  -  5  » 

C  FRECGH  IS  THE  INITIAL  FREQUENCY  (GIGAHERTZ) 

C  FIFTH  CARD-  IN  FCRMATU2) 

C  IN-INDEX  NUMBER  OF  RELATIVE  DIELECTRIC 

C  PERMITTIVITI tS  TO  BE  READ  IN  NEXT  CARDS 

C  SIXTH  AND  FOLLOWING  CARDS  CONTAIN   VALUES  OF 

C  EPR2  IN  REGION  TWO  ONE  FOR  EACH   CARD  (F10.5) 

C 

C  THE  METHOD  USED  IN  FINDING  THE  EFFECTIVE 

C  WAVELENGTH  IN  THE  STRUCTURE  IS  BASED  ON  A  RCOT 

C  SEARCHING  ALGORITHM,   SPECIFICALLY  DEVELOPED 

C  FPR   THIS   STUDY 

C  PROGRAM  DEVELOPED   BY  LT(JG)   AHMET  M.  TUFEKCIOGLU 

C  TURKISH  NAVY 

C  THESIS  ADVISOR  -PROF.  JEFFREY  B.  KNORR   PH.D. 

C  NAVAL  POSTGRADUATE   SCHOOL  MONTEREY  CALIFORNIA 

C  SEPTEMBER   1974 

C  / 

C  ALL  VARIABLES   ARE  DOUBLE  PRECISION   REAL*8 


C 


IMPLICIT  REAL*8( A-H,K-Z) 

DIMENSION  FYA(80),  FYB(80),  FYC(80),  FYD(80),  TRRA(2 

DIMENSION  YX(80),  YY(80)f  VZ(80),  YC(bO),  YF(80) 

DIMENSION  TSK2),  TS2(2),  TS3(2),  TSTK) 

COMMON  /FER/  PNS, IFSW 

COMMON  /ZEYNEP/  S,IEO 

COMMON  /ADA/  D,wOVD,WW 

COMMON  /AT/  TA, T3,TC  ,TD,PI , BETAD  , TE2 , EPR2 , IFLG 

COMMON  /EMRE/  T E , T F , F A2 , F A3 , FBI , F B2 , FB9 , FC 1 , FC2, F01 

EXTERNAL  GZP,GZIM 

DATA  TS1/'DIELECTR« ,'  FERRITE'/ 

DATA  TS2/'   SINGLES'  COUPLED'/ 

DATA  TS3/'    EVEN  • ,'    ODD   '/ 

READ  (5,18)  IDF, ISC, IEO 

PE^D  (5,19)  EPR1,MR1, MR2 

READ  (5,19)  DDfWWlfSS 

REAL  (5,20)  HREQGH 

READ  (5,21)  IN 
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c 
c 


c 
c 


PI  =  3. 14159265358979D0 

PISQ  =  PI**2 

EPS  =  0.1**1 

D  =  CD*l.E-3 

S  =  SS*l.E-3 

WW  =  WWl*l.E-3 

WCVD  =  WW/D 

FPEQDD  =  FREQGH 

SETTING  LOWER  AND  UPPER 

BCB  =  -2.*PI/WW 

BBB  =  -BCB 

DC  17  IK=1,IN 

READ  (5,20)  EPR2 

KAK  =  DSQRT(EPR2) 

FFECGH  =  FREQDD 

WRITE  (6,22) 

WRITE  (6,25) 

IF  (ISC.EQ.l) 

WRITE  (6,26) 

GC  TC  2 

WRITE  (6,23) 

WRITE  (6,24) 

IFRS  =  IDINT( 

MAXIMUM  FREQUENCY   RANGE 

IN   GIGAHERTZ   TG  AVOID 

IF  (IFRS.C-T.50)  IFRS=50 


LIMITS  OF  INTEGRATION 


TS2(ISC)  ,TS1(IDF) 
WOVD 

GO  TO  1 
WW1,SS,TS3(IE0) 

WW  1 

EPR2,MR2,DD 

.15/<KAK*D)J 

TO  BE 
HIGHER 


COVERED 
ORDER  MODES 


DC  14  IS=1,IFRS 
C      IFLG   IS  A  FLAG  WHEN  EQUALS  TO  1   CHARACTERISTIC 
C      IMPEDANCE   CALCULATED.  IFSW  IS  THE  FERRITE   SWITCH 
C      WHEN  IS  EQUAL  TU  1  FERRITE  SUBSTRATE   RELATED 
C      CALCULATIONS   ARE  DONE 

IFLG  =  0 

IFSW  =  0 

FREG  =  FREQGH*1.E9 

LAM  =  3.E8/FREQ 
C 

DCVL  =  D/LAM 
C      RCOT  FINDING  ROUTINE  STARTS 
C      USTL^UPPER  LIMIT  FCR   L/LPR 
C      ALTL=LQW£k  LIMIT   FCR   L/LPR 

USTL  =  KAK-.01 

IF  (IS- 2)  3,4,4 

3  ALTL  =  (USTL+l.J/2. 
GO  TO  5 

4  ALTL  =  LCVLPR 

5  LCVLPR  =  USTL 
ALTL1  =  ALTL 

I A  =  1 
C      INTERMEDIATE   ALGEERAIC   STEPS 

6  BETAD  =  2.0*PI*D0VL*L0VLPR 
TA1  =  240.0*PISQ*D0VL 

TA  =  TA1*MR2 
TBI  =  DOVL/60.0 
TE  =  TB1*EPR2 
TC  =  TA1*MR1 
TD  =  TB1*EPR1 
C 

DCVLSQ  =  D0VL**2 
TE1  =  LGVLPR**2 
TE2  =  4.0*PISQ*D0VLSQ 
TE  =  -TE2*(TE1-1.0) 
TF  =  TE2*( EPR2-TE1) 
FA2  "    TE/TC 
FA3  =•  TF/TA 
FBI  =  BETAD/TC 
FB2  =  BETAD/TA 
•  F69  =  TE/TF 
FC1  =  BETAC/TE 
FC1  =  TD/TE 
FC2  =  TB/TF 
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C      VALUES  OF  ALFA   WHERE   GAMMA  2   SWITCHES 
C      REAL  TO  IMAGINARY 

ALF  =  DSQRT(TF/D**2) 

ALFN  =  -ALF 
C      EVALUATION  OF   INTEGRAL 

CALL  DGG24  ( BOB , ALFN , GZR, Yl ) 

CALL  DQG24  I ALFN , ALF, GZ I M, Y2) 

CALL  DCG24  ( ALF , B6B, GZR, Y3  ) 

GZ  =  Y1+Y2+Y3 

IF  (DABS(GZ).LT.EPS)  GO  TO  11 

IF  (IA.NE.l)  GO  TO  7 

IA  =  0 

LOVLPR  =  ALTL 

GZ1  =  GZ 

GC  TO  6 

7  IF  (DABS(GZl-GZ) . GT .DABS ( GZ  1 )  )  GO  TO  9 
IF  (ALTL.NE.ALTL1)  GO  TO  8 

ALTL  =  ALTL+O.OOl 
GC  TO  10 

8  USTL  =  ALTL 

ALTL  =  (ALTL1  +  USTD/2. 
LCVLPR  =  ALTL 
GZ1  =  GZ 
GC  TO  6 

9  ALTL1  =  ALTL 

10  ALTL  =  (USTL  +  ALTD/2. 
LCVLPR  =  ALTL 

GO  TO  6 

11  YX(IS)  =  DOVL 
YY(IS)  =  1. /LOVLPR 
YZ(IS)  =  LCVLPR 
YF(IS)  =  FPEQGH 
IFLG  =  1 

C      EVALUATION  OF  INTEGRAL   FOR  CHARACTERISTIC  IMPEDANCE 

r/ui  nor.9i  'BCBjA'  FNj^g-C} 

CALL  DCG24  I  ALFN, ALF, GZIM, C2J 

CALL  0GG24  ( ALF, BBB, GZR, C3 J 

CHIMP  =  (C1+C2+C3)*(1./(8.*PI) )*D*2.0 

IF  USC.EQ.2)  CHIMP-CHIMP/2.D0 

YC(IS)  =  CHIMP 

IF  (IDF.NE.2)  GO  TO  13 
C      FERRITE  CALCULATIONS  STARTS. 
C      SET  OF  PARAMETERS   USED   FOR   SPECIFIC  CASE. 

HZERO  =  6.0D0 

WZERO  =  2.8D0*HZER0 

WZERSQ  =  WZER0**2 

WALFA    =    O.1O5D0 

WALFSQ    =    WALFA**2 

WM    =    3.36D0 

WMSQ    =    WM**2 

FtRWSQ    =    (2.*PI*FREQGH)**2 

FERDtN    =    (WZERSQ-FERWSQ-WALFSQ)**2+4.*WZERSC*WALFSQ 

XKIPRI    =     { WM*WZERO*(WZERSQ-FERWSC+WALFSQ) J/FERDEN 

XKIDPR    =    UM*WALFA*(WZERSQ+FERWSG+wALFSQ))/FERDEN 

IFSW    =     1 
C  INTEGRATION       LIMITS    DOUBLED    FOR    FERRITE     PART. 

BAB    =    2.*6CB 

PllP  =  — BAR 
C      PNS   IS  POSITIVE  NEGATIVE  SWITCH   FOR  +Z  ANC  -Z 
C      DIRECTION  TRAVELING  FIELDS. 

PNS  =  -l.ODO 
C 

DC  12  IXX=1,2 
C      EVALUATION  OF   INTEGRAL 

CALL  DQG24  ( BAB, ALFN, GZR, FER1 ) 

CALL  DGG24  ( ALFN , ALF, GZIM, FER2 ) 

CALL  DGG24  ( ALF , BUb,GZR,FER3 J 

TRPA(IXX)  =  FER1+FER2+FER3 
.   PNS  =  PNS**2 

12  CONTINUE 
C 

MRMO  =  4.D-7*PI 
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13 
14 


SABIT    =     (FREQ*D*MRM0)/(4.*CHIMP) 

FYA(IS)     =    DA6S(XKIPRI*TRRA( 1)*SABIT) 

FYB(IS)     =    DAbS(XKIPRI*TRRA(2)*SABIT) 

FYC(IS)     =    DABS(XKIDPR*TRRA(1)*SABIT) 

FYDUS)     =    DaBS(XKIDPR*TRRA(2)*SABIT) 

BETTA  =  BETAO/D 

FYA(IS)  =  (BETTA+FYA( IS))/BETTA 

FYB(IS)  =  (BETTA+FYB( IS) J/BETTA 

FYC(IS)  =  8.b86*2.*PI*FYC( I S) / ( BETTA*FYA ( I S ) ) 

FYD(IS)  =  8.686*2.*PI*FYD( I S ) / ( BETTA*FYA ( I S ) ) 

FREGGH  =  FREQGH+1.000 

CONTINUE 


C 
C 
C 


15 


16 
17 


18 
IS 

20 
21 
22 

23' 
24 


WRITE  (6,27) 
WRITE  (6,28) 

DC  15  ITR=1,IFRS 

WRITE  (6,29)  YF(ITR),YX(ITR),YZ(ITR) ,YY( ITR),YC(ITR) 

IF  (IDF.NE.2)  GO  TO  17 
WRITE  (6,30) 

DC  16  IH=1,IFRS 

WRITE  (6,31)  YF(  IH),FYA(IH) , FYB ( IH) , FYC ( IH)  ,FYD( IH) 

CONTINUE 

STOP 


26 

27 

28 
29 

30 

31 


FCRMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 

•  r  T  n  t  n 

Oirvjr 

FCRMAT 
FCRMAT 

•THICK 
FORMAT 
•W/D='  , 
!   lH*/2 
FORMAT 

56X,»S 
FCRMAT 
'L/LPR' 
FCRMAT 
FCRMAT 

F6.2,l 
FCRMAT 

•BF+/8 
FCRMAT 
END 


(311) 
(3F10.5) 

(F10.5) 

(  12) 

(• 1' ,//55X, 'MICROSTRIP  ANALYS I S » //5  5X, A8 ,2X , 

■  /  /  ^  =  Y. .  £  ?. ,  2  X  ,  '  S 0 3  ^ T  R  *\T  p  •  '  » 

'WIDTH  OF  STRIP=»  ,F5.2) 

• EPR2=' ,F5.2,2X, «MR2=' ,F5.2//56X, 

F4.1y2Xf *MM'/J 

13(1H*)/2(61X,1H*,11X,1H*/)  ,61X,1H*, 


(  '0' 
(  40» 
NESS 
(«  0' 
F4.1 
(61X 
CO' 
EPER 
CO' 
,6X, 
(•  • 
CO' 
X,  '0 
(•  0' 
PR"  , 
CO' 


,55X, 

,55X, 

D=«, 

,60X, 

»3XI 

,  IH*  , 

,55X, 
AT  I  ON 
,32X, 
•LPR/ 
i30X, 
,32X, 
HM  »  /  ) 
,8X,« 
10X,  • 
,3X,F 


11X,1H*/) ,61X,13(1H*)///) 

•WIDTH    OF    STRIP=« ,F4.1,2X, 
• ,F4. 1,2X, «MM«//55X,A6,2X, 
•FREQUENCY     ',2X,«     D/LAMBDA     », 
L« ,5X, «CH    IMPEDANCE' ) 
5(2X, 10(1H_) )//) 
Fo.o,lX,'GHZ',3(2X,F10.5),2X, 


• NM'// 

•  NODE' 


5X 


FREQUENCY' ,8X, • BF-/BPR  ',10X, 
2*PI*8.668*DELALFA/B  **NEG  AND  PCS') 
10.3,4(4X,F15.6) ) 
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REAL    FUNCTI0NGZR*8(ALFA) 

ROUTINE    DEFINES    THE       POLYNOMIAL    FOR    HYPEBOLIC    CASE 

IMPLICIT    REAL*8( A-H,K-Z) 

COMMON  /FER/  PNS,  IFSW 

COMMON  /ADA/  D,WCVD,WW 

COMMON  /AT/  TA,TBTTC,TDtPI , BETAD , TE2 , EPR2 » I FLG 

COMMON  /EMRE/  T E , T F  , FA2 , FA3 , FB 1 , FB2 , FB9 , FC 1 ,FC2» FD1 

ALFAD  =  ALFA*D 

TGI  =  ALFA0**2 

TG  =  TG1-TE 

VG  =  DSQRT(TG) 

B6  =  FB1*ALFAD/VG 

EC  =  FC1*ALFAD 

TH  =  TG1-TF 

VH  =  DSQRT(TH) 

FA1  =  l./DTANHCVH) 

BA  =  FA1*FB2*ALFAD/VH 

FFF1  =  -(FA2/VG+FA3*FA1/VH) 

FFF2  =  6B+BA 

FFF3  =  -FFF2 

FFF4  =  eC*(6B+8A*FB9)-(FCl*VG+FC2*FAl*VH) 

DENOM  =  FFF1*FFF4-FFF2*FFF3 

GG4  =  -FFF1/DENOM 

GG2  =  FFF2/DENOM 

CALL  TRAN  (  ALFA , ALFAD, GX) 

IF  (IFLG.EC.l)  GO  TO  1 

GZR  =  GG4*GX 

RETURN 

RCAA  =  GG4/TE 

RCAA1  r  ALFA0*FD1*GG4 

RCAB  =  <1./(TC*VG) )*(GG2-RCAA1) 

PRA  =  GX*RCAA**2 

KRB  =  GX*RCAB**2 

RCA    =    RCAA*RCA6*GX 

RCAC    =    GG4/(TF*DSINH( VH) ) 

RCAD    =     l./(TA*VH*DSIi\H{VHH 

RCAE  =  -GG2+RCAA1*FB9 

RCB  =  RCAC*RCAD*RCAE*GX 

RRC  =  GX*RCAC**2 

RRD  =  GX*<RCAD*RCAE)**2 

P1B  =  ( ITG1+TG)/VG)*(TD*RRA+TC*RRB)*BETAD 

PIC  =  2.*ALFAD*(BETAD**2+TE2)*RCA 

PUNT  =  Pie+PIC 

RKA1  =  DSINH(2.*VH) 

RKA  =  RKA1-<2.*VH) 

RKB  =  TG1*BETAD*TB*RRC/VH 

RKC  =  BETAD*TA*VH*RRD 

TAVTB  =  TE2*EPR2 

RKD  =  ALFAD*((TAVTB)+BETAD**2)*RCB 

RKE    =     RKB+RKC-RKD 

RKG    =    RKA1+(2.*VHJ 

RKH    =    TG1*BETAD*TA*PRD/VH 

RKI    =    BETAD*TB*VH*RRC 

RKK    =     ITAVTB+BETA0**2)*ALFAC*RCB 

P2INT    =    RKA*RKE+RKG*(RKH+RKI-RKK) 

IF     (IFSW.EQ.l)     GO    TO    2 

GZR    =    P1INT+P2INT*0.5 

RETURN 

FERAR  =  RKA/(2.*VH) 

FERB  =  {-TG1*BETAD**2+TF**2)*RRD 

FERC  =  PNS*T6*BEIAD*ALFAD*VH*2.0*RCB 

FERD  =  TH*RRC*TB**2 

GZR  =  FERAR*(FERB+FERC-FERD) 

RETURN 

END 
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DS 
FB 
FD 
-T 
DS 

:   D 

FA 


REAL  FU 
IMPLICI 
CCMMON 
COMMON 
CCMMGN 
COMMON 
ALFAO  = 
TGI  =  A 
7G  =  TG 
VG  = 
BB  = 
BC  = 
TF  = 
VH  = 
FA1 
BA  = 
FFF1  = 
FFF2  = 
FFF3  = 
FFF4  = 
DENOM  = 
GG4  =  - 
GC-2  =  F 
CALL  TR 
IF  (IFL 
GZIM  = 
RETURN 
RCAA  = 
RCAA1  = 
RCAB  = 
RPA  =  G 
RRB  =  G 
RCA  =  R 
RCAC  = 
RCAD  = 
RCAE  = 


NCTI0NGZIM*8(ALFA) 

T  REAL*8( A-H,K-Z) 

/PER/  PNSi IFSW 

/AT/  TA,TB,TC,TD,PI , BETAO  ,TE2 , EPR2 , 1 FLG 

/ADA/  D,WCVD,WW 

/EMRE/  TE,TF,FA2,FA3,FB1,FB2,F69,FC1,FC2,FD1 

ALFA*D 
LFAD**2 
1-TE 
QRT(TG) 
1*ALFAD/VG 
1*ALFAD 
Gl+TF 
QRT(TH) 
COTAN(VH) 
i*FB2*ALFAD/VH 
-(FA2/VG-FA3*FA1/VH) 
BB-6A 
-FFF2 
BC*(3B-BA*FB9)-(FC1*VG+FC2*FA1*VH) 

FFF1*FFF4-FFF2*FFF3 
FFF1/DENGM 
FF2/OENOM 

AN  (ALFA, ALFAO, GX) 
G.EC.l)  GO  TO  1 
GG 


4*GX 


RCB 
RRC  = 
PRO  = 
P1B  = 
PIC  = 
PUNT 
RKA1  = 
RKA  = 
RKB  = 
RKC  = 
TAVTB 
RKD  = 
RKE  = 
RKG  = 
RKH  = 
RKI  = 
RKK  = 
P2INT 
IF  (IPS 
GZIM  = 
RETURN 
FERAI  = 
FERB  = 
FERC  = 
FERD  = 
GZIM  = 
RETURN 
END 


GG4/T 

ALFA 
(l./( 
X*RCA 
X*RCA 
CAA*R 
GG4/( 
I .  /  ( T 
GG2-R 
CAC*R 
X*RCA 
X*(RC 
(TG1  + 
.*ALF 

PI  6+ 
DSINt 
.*VH- 
G1*BE 
ETAD* 

TE2* 
LFAD* 
KB  +  RK 
KAl  +  ( 
G1*BE 
ETAD* 
TAVTB 

RKA* 
w.  EG. 
PUNT 


E 

D* 

TC 

A* 

B* 

CA 

TF 

A* 

CA 

CA 

C* 

AD 

TG 

AD 

PI 

2. 

RK 

TA 

TA 

EP 

(  ( 

C+ 

2. 

TA 

TB 

+  B 

RK 

1) 

+  P 


FD 
*V 
*2 
*2 
B* 
*D 
VH 
Ai 
D* 
*2 
*R 
)/ 
*( 

C 

*v 

Al 

D* 
*V 
R2 
TA 
RK 
*V 
D* 
*V 

ET 
E+ 

21 


1*GG4 

G)  )*(GG2-RCAA1) 


C  X 

SIN(VH) ) 
*  n  <;  t  m  t  \j  h  \  ) 

*F69 
RCAE*GX 


r  a  f  i  #  #  ? 
vg)*Ttd*rra+tc*rrbj*betao 

betad**2+te2)*rca 


H) 


TB*RRC/VH 
H*RRD 


VTB)+BETAD**2)*RCB 

D 

H) 

TA*I 

H*Ri 

AD: 

RKG: 

0    Tl 

NT*' 


;RRD/VH 

;RC 

=  *2)*ALFAD*RCB 

l*( PKH+RKI-RKK) 

0    2 

0.5 


(RKAI+2.*VH)/(2.*VH) 
(-TG1*BETAD**2+TF**2)*RRD 
PNS*TB*3ETA0*ALFAD*Vh*2.0*RCB 
TH*RRC*T5**2 
FERAI*(FER6+FERC-FERD) 
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SUBROUTINE  TRAN  ( ALTR , ALFD, GXTR  ) 

IMPLICIT  REAI_*8(A-H,K-Z) 

CCMMON  /ZEYNEP/  S , I EO 

COMMON  /ADA/  D,WOVD,WW 

CB  =  ALFD*hCVD*0.5 

GX1  =  ( (l./CB)*(DSIN(CB)))**2 

IP  (  IEO.EQ.O)  GO  TO  2 

CD  =  ALTR*(S+WW)*0.5 

IF  (IEQ.EQ.21  GO  TO  1 

GX1  =  GX1*(2.*DC0S<CD) )**2 

GC  TO  2 

1  GX1  =  GX1*(2.*DSIN(C0) )**2 

2  GXTR  =  GX1 
RETURN 
END 
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